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DIFFERENTIAL 
EQUATIONS 


EQUATIONS WITH SEPARABLE VARIABLES 


A first order differential equation is an equation involving an independent variable f, 
a dependent variable y, and the derivative dy/dt. In many applications, the independent 
variable r is time. A first order differential equation can be put in the following 
form, where f(t, y) is continuous in both ¢ and y. 


FIRST ORDER DIFFERENTIAL EQUATION 


dy 
(1) pate) 


For instance, if f is time and y is the position of a particle at time f, the differential 
equation (1) gives the velocity of the particle in terms of time and position. A 
differential equation gives information about an unknown function y(t). The general 
solution of a first order differential equation is the family of all functions y(t) that 
satisfy the equation. Each function in this family is called a particular solution of 
the differential equation. In most cases, the family of functions will depend in some 
way ona constant C, and the graphs of these functions will form a family of curves 
that fill up the (t, y) plane but do not touch each other, as in Figure 14.1.1. 

Some examples of first order differential equations were solved in Section 8.6. 
For instance, it was shown that the general solution of the differential equation 


dy 
—_—= y(1—-—y 
"a 0 ae) 
: 1 
Is yt) = Ta Cer y(t) = 0. 


There is one particular solution for each value of the constant C, and one additional 
particular solution )(t) = 0. The graph of this general solution is shown in Figure 
14.1.2. 
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Figure 14.1.1 


Figure 14.1.2 


In most applications, a first order differential equation will describe a 
process that starts at some initial time tg. In order to determine a particular solution, 
we need both the differential equation and the value of y(t) at the initial time to. A 
first order initial value problem is a pair of equations consisting of a first order 
differential equation and an initial value. 
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FIRST ORDER INITIAL VALUE PROBLEM 


(2) a =f(t.¥), — Wto) = Yo- 
t 

An initial value problem usually has just one solution, which will be a particular 
solution of the differential equation. This can be seen intuitively as follows. Figure 
14.1.3 shows a moving point controlled by an infinitesimal driver with a steering 
wheel. At tf = to, the moving point starts at (to, yo). At each ¢ > fo, the infinitesimal 
driver measures his position (f, y), computes the value of f(t, y), and turns the steering 
wheel so that the slope will be f(t, y). The curve traced out by this point will be the 
solution of the initial value problem. 


ty Yo) 


Figure 14.1.3 


The general problem of solving a differential equation is difficult. In Section 
8.6, we gave a method for solving differential equations of an especially simple 
type, called equations with separable variables. As a starting point in this chapter, 
we take another look at these differential equations. 


FIRST ORDER DIFFERENTIAL EQUATION WITH SEPARABLE VARIABLES 


dy 

(3) To g(t)h(y). 

In an equation with separable variables, dy/dt is a product of a function of t and a 
function of y. In particular, dy/dt = g(t) is a differential equation with separable 
variables in which h(y) is the constant 1. Similarly, dy/dt = h(y) is a differential 
equation with separable variables in which g(t) is the constant 1. An equation with 
separable variables can be solved by separating the variables and integrating both 
sides of the equation. 


METHOD FOR SOLVING A DIFFERENTIAL EQUATION WITH SEPARABLE VARIABLES 


dy _ 
ae g(t)h(y). 


14.1 EQUATIONS WITH SEPARABLE VARIABLES 


Step 1 Find all points y, where h(y,) = 0. For each such point, the constant function 
y(t) = y, is a particular solution. 


Step 2. Separate the variables by dividing by h(y) and multiplying by dt, 


ra dy = g(t) dt, 
then integrate both sides of the equation. That is, find antiderivatives of each 
side. 
KW)=p5, GO=E0, 
Ay) 
so that K(Qy) = GI) + C. 


If possible, solve for y as a function of t. 


Step 3 The general solution is the family of all solutions found in Steps 1 and 2. It 
will usually depend on a constant C. 


Step 4 If an initial value y(to) = yo is given, use it to find the constant C and the 
particular solution of the initial value problem. 


Remark The cases h(y) = 0 and h(y) # 0 must be done separately in Steps 1 and 2, 
because the division by h(y) in Step 2 cannot be done when h(y) = 0. 


The general solution of a differential equation dy/dt = g(t), where dy/dt 
is a function of f alone, is just the indefinite integral 


jee feo dt = G(t) + C. 


In this case, C is the familiar constant of integration, which is added to a particular 
solution. For example, the general solution of the differential equation dy/dt = 1/t 
isy =In|t|+C. 

In the examples that follow, the constant C appears in a more complicated 
manner. 


EXAMPLE 1 Solve the initial value problem 
dy _ 
dt 

Step 1 —2y =0 when y = 0. Thus the constant y(t) = 0 is a particular solution. 


“2y5. Vl); 


Step 2 Separate the variables and integrate both sides. 


la —2 dt, 
y 

In|y| = —2¢ + B, 
ly| = eFe7 74, 
y= Ce, 


where C = e? ify > 0,andC = —e® ify < 0. 


849 


850 


Step 3 


Step 4 
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General solution: 
y(t) = Ceo, 
where C is any constant (C can be 0 from Step 1). 


Substitute 1 for t and —S5 for y, and solve for C. 


—5 = Ce™*", —5=Ce?, C = —Se?, 


Particular solution: 
y(t) = —Se?e~ 2 = — Sel? 2”, 


The graph of the solution is shown in Figure 14.1.4. 


(1,~5) 


Figure 14.1.4 Example 1 


EXAMPLE 2 Find the general solution of the differential equation 


Step 7 
Step 2 


Step 3 


Step 4 


“ = y’ sine. 
Then find the particular solution with y(0) = 4. 
y? = 0 when y = 0. Thus y(t) = 0 is a constant solution. 
y ?dy =sintdt, 
—y '= —cost+C, 
y=(cost — C)7t. 
General solution: 
y(t) = 0 and y(t) = (cost — C)7'. 
0) =3 =(cos0- Cyt = (1-0, 
Po Ee OS, 


~ 


Particular solution: 
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y(t) = (cost + 1)7h. 


The particular solution to Example 2 is illustrated in Figure 14.1.5. It is 
defined only for —x < t < x and approaches oo as t approaches 7. It is said to have 
an explosion at t = 7. 


Figure 14.1.5 Example 2 


To avoid errors, the general solution can be checked by differentiating. 
The solution for Example 2 is checked as follows. 


y(t) = (cost — C)71. 


as required. 


a —(cost ~ C) ?(—sin t) 


dt 


= (cost — C)~ (sin t) 


= y’sint, 


PROBLEMS FOR SECTION 14.1 


In Problems 1-12, find the general solution of the given differential equation. 


1 
3 
5 


7 


9 
11 


In Problems 13-18, solve the initial value problem. 
y(1) = 2 14 


13 


y’ =t-sin(t?) 


y =e” 

y= y? -— 1 
yoy 
y=(1t+y’)e' 
y =ytant 


y = ye, 


2 
4 
6 


8 


10 
12 


—3t 


y=e 
y=y 

y =tyy + 1) 
_ ti 
[ = Oya A 


y =./1 — y* cost 


y =ysint 


y=n/y, 


y(0) = 3 
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n 
15 y=—, y(1) = -2 16 yo=ty—y+2r—2, y(0) = 0 


¥ 
He y=Q7-3y+2)4, yU)=2 
18 y=e(y — 4)sine, y(2) = 4 


FIRST ORDER HOMOGENEOUS LINEAR EQUATIONS 


In this section we study the following special type of differential equation. 


FIRST ORDER HOMOGENEOUS LINEAR DIFFERENTIAL EQUATION 


(1) y’ + p(t)y = 0. 
It is understood that ¢ varies over some interval in the real line, and p(t) is a con- 
tinuous function of t in the interval. The equation is called linear because y and y’ 
occur only linearly and homogeneous because the right side of the equation is zero. 
The equation 

yoky, or y—ky=0, 


for exponential growth (see Section 8.6) is an example. The first order homogeneous 
linear differential equation (1) has separable variables, because it can be written as 


dy 
ae — p(t)y. 


Its solution is given by the next formula. 
METHOD FOR SOLVING FIRST ORDER HOMOGENEOUS LINEAR DIFFERENTIAL 
EQUATION (1) 


The general solution is 


y(t) = Ce, 
where P(t) is an antiderivative of p(t). That is, 
y(t) = Ce I POa, 


This formula is obtained by the procedure described in Section 14.1 for differential 
equations with separable variables, as follows. First write the equation in the form 


dy 
7 Ties — p(tyy. 


Step 7 There is a constant solution y(t) = 0. 
Step 2 Separate the variables and integrate: 


yp! dy = —p(t) dt. 


In|y| = - | pear +B 
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Now solve for y. 
—Jf pit) dt+B 
» 


lyl =e 
y= Ce 7S Pa)ar 


where C =e? if y>0,and C = —e® ify <0. 


Step 3 Combining Steps 1 and 2, we get the general solution 


y(t) = Cem! Poa, 


Remark The case C = 0 gives the constant solution y(t) = 0 of Step 1. 


Discussion The constant of integration in the indefinite integral 


f ve at 


will be absorbed in the constant C. 

The particular solution for the initial value y(t9) = yo is found by substituting 
and computing C. Notice that any two particular solutions of the same 
homogeneous linear differential equation differ only by a constant factor. 
If x(t) is any nonzero particular solution, then the general solution is Cx(t). 


EXAMPLE 1 


(a) Find the general solution of the equation y’ + ycost = 0. 
(b) Find the particular solution with initial value y(0) = 3. 


(c) Find the particular solution with initial value y(2) = 3. 


SOLUTION 
(a) First evaluate the integral 


{ cos ae = sint + B. 

General solution: 

y@) = Ce", 
(b) First substitute and solve for C. 

y(0) =$ = Ce 9 = CL =C. 

Particular solution: 

(0) = Feosn|, 
(c) Substitute and solve for C. 


y(2) = 
Cc 


— —sin 2 
= Ce on 


2 
Losin? _ 1.2413, 


I 


Particular solution: 
y(t) = 1.241327 si8 


The solution to this example is shown in Figure 14.2.1. 
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(0, 4) 


Figure 14.2.1. Example 1 


EXAMPLE 2 
(a) Find the general solution of the equation ty’ + 3y = 0 fort > 0. 
(b) Find the particular solution with the initial value y(1) = 2. 


SOLUTION 
(a) We first put the equation into the homogeneous linear form (1) by 
dividing by t: 


y + 3t7'y =0. 
Next evaluate the integral, 
facta 3Int + B. 
The constant of integration B is absorbed into the constant C, and the 


general solution is 
ij = te aC, 


(b) The particular solution with initial value y(1) = 2 is 
y(t) = 2173. 
The solution to this example is shown in Figure 14.2.2. 


Ay 


(1,2) 


1 


Figure 14.2.2 Example 2 


14.2 FIRST ORDER HOMOGENEOUS LINEAR EQUATIONS 


The coefficient 3t~ 1 in the equation 
y + 3t7'y =0 


of Example 2 is discontinuous at t = 0. However, it is continuous on the interval 
t > 0 and on the interval t < 0. An initial value at a positive time t) > 0 will 
determine a particular solution only for the interval t > 0, while an initial value 
at a negative time tg < 0 will determine a particular solution for the interval t < 0. 
Each interval must be solved separately. The next example is like Example 2 but 
on the negative time interval. 


EXAMPLE 3 


(a) Find the general solution of the equation y’ + 3t~'y = 0 from Example 
2, but for the interval t < 0 instead of t > 0. 


(b) Find the particular solution of the initial value problem with y(—2) = 1. 


SOLUTION 
(a) This time we integrate with a negative f, 


[setae=3in|e} +B =3In(-) 4B. 


The general solution for t < 0 is thus 
y(t) = Ce F 8-9 = C(—-F) 3 = — Cr 3, 
or y(t) = At” 3, 


where A is the constant —C. 


(b) The particular solution with the initial value y(—2) = 1 is found by 
solving for A: 


1=A(-2)°37, A=-8, yt) = —873. 


The solution to this example is shown in Figure 14.2.3. 


Figure 14.2.3 Example 3 


For some purposes, it is useful to describe the solution of a differential 
equation using a definite integral from some point a to t, instead of using an indefinite 
integral. In the definite integral form, the general solution of the linear homogeneous 
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differential equation 


. y+ pOy = 0 
: y(t) _ Ce 7's p(s) -_ 
where a is any point in the interval I, and C is any real number. This formula is helpful 


in a problem where one cannot evaluate the integral of p(r) exactly and must use a 
numerical approximation. The formula holds because the integral 


[r) ds 


is an antiderivative of p(t) by the Fundamental Theorem of Calculus. The choice 
of the endpoint a does not matter because a change in the value of a will be absorbed 
by a change in the value of the constant C. If we are given an initial value y(to) = yo, 
the particular solution can again be found by substituting and solving for C. If we 
take a = fo, the constant C will be equal to yg, 


Yo = y(a) = Ce $4 PO) ds — Ce® = C, 
Thus the particular solution of the initial value problem (1) with y(a) = J is given by 


y(t) = yoe 4 p(s) ds’ 


PROBLEMS FOR SECTION 14.2 


In Problems 1-4, find the general solution of the given differential equation. 
1 y+ 5y =0 2 yi 2y=0 
3 y+, =0 4 try =0 


In Problems 5-12, find the particular solution of the initial value problem. 


5 yry=0 yQ=4 6 y-3y=0, yl) = -2 
7 y + y'sint = 0, y(t) = 1 8 y + ye =0, yO) =e 
9 rtyJ/l+=0,  y(0)=0 
10 y” + ycos (e') = 0, 3(0) = 0 
11 ty’ — 2 =0, yl) = 4, t>0 
12 Pyty=0 yD=-2 +60 
13 ry’ = 2y, yQ) = 1, t>0 
14 fy =2y, w=0, 1>0 
15 A function }(t) is a solution of the differential equation 
yo tky =0 
for some constant k. Given that (0) = 100, and )(2) = 4, find k and find y as a function 
of t. 
16 A function y(t) is a solution of the differential equation 
y tty =0 
for some constant k. Given that )(0) = I,and )(1) = e7! 3, find k and find yas a function 
of t. 
17 A bacterial culture grows at a rate proportional to its population. If it has a population 


of one million at time ¢ = 0 hours and 1.5 million at time ¢ = 1 hour, find its population 
as a function of r. 
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18 A radioactive element decays with a half-life of 6 years. Starting with 10 Ib of the element 
at time t = 0 years, find the amount of the element as a function of f. 


FIRST ORDER LINEAR EQUATIONS 


We shall now give a method for solving a differential equation of the following type. 


FIRST ORDER LINEAR DIFFERENTIAL EQUATION 


(1) y + pOy=f(. 


Both p(t) and f(t) are continuous functions of t, where t varies over some interval J 
in the real line. When f(t) is the constant function with value 0, the equation is a 
homogeneous linear differential equation of the type studied 1n Section 14.2. 

First order linear differential equations arise in models of population 
growth with immigration. Suppose a population y(t) has a net birthrate of b(t) and 
net immigration rate of f(t). The net birthrate b(t) is the excess of births over deaths 
per unit of population in one unit of time. In a small period of time of length Ar, 
the difference of births and deaths is b(t) + y(t)+ At, and the net immigration is 
f(t)+At. Then the population will be a solution of the differential equation 
y = b(t)y + f(t), which is the same as equation (1) with p(t) = — b(t). 

The size of a bank account that earns interest and also changes due to 
deposits and withdrawals can be described by a first order linear differential equation. 
If the account earns interest at the rate of r(t) at time ¢, and the net deposit per unit 
of time is f(t), then the account size y(t) will be a solution of the differential equation 
(1) with p(t) = —r(t). 

The next theorem will be helpful in solving an equation of the type (1). 


THEOREM 1 
Suppose that y(t) is a particular solution of the first order linear differential 
equation 

(1) y+ poy =f, 
and x(t) is a nonzero particular solution of the corresponding homogeneous 
equation 

(2) x’ + p()x = 0. 


Then the general solution of the original equation (1) is 


y(t) + Cx(e). 


We already know from Section 14.2 how to solve the homogeneous linear 
equation (2). So if we can find one particular solution of the linear equation (1), we 
can use Theorem 1 to find the general solution. We postpone the proof of Theorem I 
to the end of this section. 

A particular solution of a linear equation (1) can be found by the method 
called variation of constants. Start with a particular solution x(t) of the corresponding 
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homogeneous equation (2). For any constant C, Cx(t) is also a solution of (2). Now 
replace the constant C by a variable v(t), and see what happens. Let (f) = v(t)x(t). 
We shall compute the left side of equation (1), y’ + p(t)y. If it turns out to be equal to 
f(t), then y(t) will be a particular solution of (1) as required. We carry out the com- 
putations using the Product Rule for derivatives. 


y = ox, 
y+ py = (vx)! + pox 
= v'x + vx’ + pox 
= v'x + v(x’ + px). 
Since x is a solution of the homogeneous equation (2), 
x’ + px = 0. 
Therefore y’ + py = vx. 
Thus if we can find a function v(t) such that 
v@x(t) =f, 
then y(t) = v(£)x(t) 
is a particular solution of the linear equation (1). 


Putting all the ideas together, we have a method for solving a first order 
linear differential equation. 


METHOD FOR SOLVING A FIRST ORDER LINEAR DIFFERENTIAL EQUATION 


(1) y+ py =f. 
The corresponding homogeneous linear differential equation is 
(2) x’ + p(t)x = 0. 


Step 1 Find a nonzero particular solution x(t) of the corresponding homogeneous 
linear differential equation (2). By the method of Section 14.2, we may take 


x(t) = eS PO at 


Step 2. Find a function v(t) whose derivative is given by 


v(t) = a. 
This is done by integration, 
v(t) = o dt. 


Step 3 The general solution of (1) is 
y(t) = v(t)x(t) + Cx(Z). 


Step 4 If an initial value is given, the particular solution for the initial value problem 
is found by substituting and solving for the constant C. 


Discussion Step 2 gives us a function v(t) for which v’()x(t) = f(t). Therefore, by 
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our previous discussion, v(t)x(t) is a particular solution of the linear 
equation (1). 
Step 3 is then justified by Theorem 1. 


EXAMPLE 1 __ Find the general solution of the equation 
yt+3tty=, t>0. 
Then find the particular solution with the initial value y(1) = 3. 
Step 7 The corresponding homogeneous equation is 
x’ + 3t7'x =0. 


From Example 2 in Section 14.2, a particular solution is 


x=t3. 
t? 
Step 2 v= pox) = 
1 
= {5}0°. 
: z} 
Step 3 The general solution is 
y = vx + Cx, 
or y=@0+Cr3, 
Step 4 W)=7=@Pr4+c1%, 
C= 


The required particular solution (Figure 14.3.1) is 
y= Qe + Ge 


a, 2 


Figure 14.3.1 Example 7 


EXAMPLE 2. A population has a net birthrate of 2% per year and a net annual 
immigration rate of 100,000sint. At time t = 0 years, the population is 
1,000,000. Find the population as a function of t. 
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The given verbal problem can be expressed as the initial value problem 

” = 0.02y + 100,000 sin ¢, y(0) = 1,000,000. 
We first put the equation in the usual form with all the } terms on the left 
side, 


y’ — 0.02) = 100,000 sin t, ¥(O) = 1,000,000. 


The corresponding homogeneous equation is 
x’ — 0.02x = 0. 


The particular solution is 
x(t) = 00-028 


yy — 100,000 sin ¢ 


e002! 


= 100,000 sin re~ 9-07", 


vu can now be found by integration by parts. With vu =sint and dw = 


e992 dt we have w = —50e7 °°! and 


[sin te °°! dr = —S50sin te~°-°*' + 50 { cos fe 02 gy, 
Similarly, 


i cos te °°?" dt = —S0 cos te °°? — 50 | sin te~ 9-97" dt, 


Combining the last two equations and solving for the integral of sin te~ °-°**, 


we get 

al - 0.028 . 

2501 e (50 sin t + 2500 cost] + K, 
— 100,000 

(=~ s504 


| sin te~ 9-97" dt = 


e~ °-°217'50 sin ¢ + 2500 cos £]. 


The general solution is y = vx + Cx, or 


100,000 
2501 


y(t) = — [50 sin t + 2500 cos tr] + Ce?! 


Substitute at f = 0. 


100,000 ‘ 
1,000,000 = — xO [50 sin 0 + 2500 cos 0] + Ce® 
100,000 
as {2500] + C, 
C = 1,099,960. 
The particular solution (Figure 14.3.2) is then 
yt) = — HOON (50 sin t + 2500 cost] + 1,099,960e°-°”, 
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107 |¥ 


10° 


100 


Figure 14.3.2 Example 2 


EXAMPLE 3 Find the general solution of the equation 
y'— sy = Ke", 
where r,s, and K are constants. 
Step 1 The corresponding homogeneous equation is 
x’—sx=0. 
It has the particular solution 


x(t) = e”. 


rt 


= Ken, 


Step 2 v(t) =—; 
e 


There are two cases,r # sandr=s. 


Step 3 (Case 71) +r#S. 


v(t) = | Ket? at 


v(t) = ae es 
r—s 


The general solution y = vx + Cx in this case is 


K 
y(t) = are + Ce™” 


Step 3 (Case 2) r=-s. In this case v(t) = K, and v(t) = Kt. The general solution 


in this case is 
y(t) = Kte* + Ce. 
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We now return to the general first order linear differential equation (1). 
Using definite integrals, we can get a single formula for the solution of equation 


(1) y+ poy =f 


by combining Steps | to 4. For Step 1, choose an initial point a, and get a particular 
solution of the corresponding homogeneous equation, 


x(t) = e7 Ja pls) ds 
For Step 2, write v(t) as a definite integral from a to ¢, 
" f(s) 
a X(S) 


Step 3 shows that the general solution is y = vx + Cx, and the final formula is 
found by substituting for v and x. 


v(t) = ds = [ f (sje?) & ds, 


GENERAL SOLUTION OF EQUATION (1), DEFINITE INTEGRAL FORM 


y(t) = e7 & p(r) ar i f (sel# py dr ay 4 c| . 


By taking t = a in the above equation, we find that C = y(a). Thus the particular 
solution of equation (1) with the initial condition y(a) = yo is obtained by replacing 
C by yo. This formula is useful when one or both of the integrals cannot be evaluated. 
In a simple problem, it is better to use Steps | to 4, which break the solution process 
into smaller parts. 

In the following example, we are able to evaluate the first integral but not 
the second, so the solution is left in a form with one definite integral. 


EXAMPLE 4 Find the general solution of the equation 
y + ycost=t. 


Step 71 From Example 1 in Section 14.2, the corresponding homogeneous equation 
has the particular solution 


—sint 


x=e 


Step 2. The function v(t) is expressed by an integral. 


if 
v= [serms ds. 
a 


We cannot evaluate the integral, so we leave it in this form. It does not 
matter which value is chosen for the lower endpoint in the integral, so we 
take the lower endpoint zero. 


Step 3 The general solution is 


y = vx + Cx, 
t 
or y(t) _ e sint [ser ds Hs Cer sink 
0) 
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EXAMPLE 5 Find the solution of the initial value problem 
y + ylnQ4+ cost) =, y) = 4. 


We are not able to evaluate the integral of In (2 + cost), so we shall use 
the definite integral form of the solution. With the initial point a = 1, the 
solution is 


t 
y(t) = ei In (2 +cos r) dr | se In (2+ cos r) dr ds By al 


1 


We conclude this section with a proof of Theorem 1. The proof uses the 
Principle of Superposition. 


PRINCIPLE OF SUPERPOSITION (First Order) 

Suppose x(t) and y(t) are solutions of the two first order linear differential 
equations 

x’ + p(t)x =f (0), 

¥ + pOy = go). 
Then for any constants A and B, the function 

u(t) = Ax(t) + By(t) 

is a solution of the linear differential equation 


u’ + p(t)u = Af(t) + Bg(t). 


Notice that all three differential equations have the same p(t). The Principle 
of Superposition follows from the Constant and Sum Rules for derivatives: 
u’ + p(t)u = (Ax + By)’ + p(t)(Ax + By) 
= Ax’ + By’ + Ap(t)x + Bp(t)y 
= A(x’ + p(t)x) + Bly’ + py) 
= Af(t) + Bait). 


PROOF OF THEOREM 7 Weare given that y and x are solutions of 


(D y+ pOy =fO 
and 
(2) x’ + p(t)x = 0. 
We must prove that a function u(t) is a solution of 
(3) u’ + p(tju = f(t) 


if and only if wu = y + Cx for some constant C. 
Assume first that u = y + Cx. By the Principle of Superposition, 
ui + pitu=f()+C-0=f(0), 


so u is a solution of (3). 
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Now assume that w is a solution of (3). Using the Principle of Superposition 
again, 
(u— y! + POU -— yb =fO —f(/) = 0. 


Thus uw — y is a solution of the homogeneous linear equation (2). The 
general solution of equation (2) is Cx. Therefore for some constant C, 


u—y=Cx and u=yt Cx. 


PROBLEMS FOR SECTION 14.3 


In Problems 1-10, find the general solution of the given differential equation. 


1 
3 
5 
7 
9 


10 


y +4y =8 2 y —2y =6 

y tty = St 4 y tey = —2e! 
y-y=P 6 2rtyst 

ty’ — 2p = 1/t, t>0o 8 ty ty= vo t>0 
ycosf+ysint = 1, —nf2<t < 1/2 

y’ + ysect = tant, —nj2<t<n/2 


In Problems 11-14, find the general solution using the definite integral form when the integral 
cannot be evaluated. 


11 
13 
15 


16 
17 


18 


19 


y+ ysint=t 12 y’ + yt? = tant, —nj2<t<n/2 
y’ + ycos(e') = 1 14 y tyel"=2e, t>0 
A population has a net birthrate of 2.5% per year and a net annual immigration equal 


to 10,000r — 40,000, where f is measured in years. At time t = 0, the population is 
y(O) = 100,000. Find the population as a function of t. 


Work Problem 15 if the net annual immigration is 1,000(cos t — 1). 


A bank account earns interest at the rate of 10% per year, and money is deposited con- 
tinuously into the account at the rate of 5? dollars per year. The earnings due to interest 
are also left in the account. If the account had $5000 at time t = 0 years, find the amount 
in the account at time t = 10 years. 


Work Problem 17 if there are no deposits but money is withdrawn continuously from 
the account at the rate of S¢? dollars per year. 


Use differential equations to prove the capital accumulation formula in Section 8.4. 
The formula says that if money is deposited continuously in an account at the rate 
of f(t) dollars per year, and the account earns interest at the annual rate r, and there 
are zero dollars in the account at time t = a, then the value of the account at time 
t = b will be 


b 
yb) = f foe ae 


EXISTENCE AND APPROXIMATION OF SOLUTIONS 


This section deals with arbitrary first order differential equations. It is optional and 
therefore can be omitted if desired. Most first order differential equations cannot 
be solved explicitly. However, it is possible to approximate a solution by a method 
similar to the Riemann sum for the definite integral. The Euler approximation 


14.4 EXISTENCE AND APPROXIMATION OF SOLUTIONS 


starts by dividing the interval [a, co) into small subintervals of length At. When 
At is real, it gives an approximate solution that can be computed numerically. When 
At is infinitesimal, it leads to a precise solution and is useful because it shows that 


a solution exists. 
Throughout this section, we shall work with a first order differential equation 


with an initial value 
(1) y=fty, ywa=yo. 


We assume once and for all that f(t, y) is continuous for all t and y. 


DEFINITION 


Let At be positive, and partition the interval [a, 00) into subintervals of 
length At. The Euler approximation for the initial value problem (1) is the 
function Y(t), a < t, defined as follows. Start the graph of Y(t) at the point 
(a, Vo). Then move from (a, yo) to (a + At, Y(a + At)) along a straight line 
with slope f(a, Y(a)). Once the value Y(t) is computed for a partition point 
t=a+kAt, move from (t, Y(d) to the next partition point (t + At, Y(t + Ad) 
along a straight line with slope f(t, Y(2)). 


The graph of Y(t) is the broken line shown in Figure 14.4.1. Each piece has 
the slope required by the differential equation (1) at the beginning of the subinterval. 
If At is small, then since f(t, y) is continuous, the slope of Y(t) should be close to 
the correct slope. Thus we would expect Y(t) to be close to a solution of (1). 


Ay 


At 


Figure 14.4.1 


The values of Y(t) at the partition points can be computed by an iteration 
that can easily be carried out on a computer. The first three values are 
Y@ = Yo, 
Y(a + At) = Yo + f(@, Yo) At, 
Y(a+2At)= Y(a+ At) + f(a + At, ¥Y(a + Ad) At. 
Given the value Y(t) for a partition point t = a + k At, the next value Y(t + At) is 
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given by the rule 
Y(t + At) = YC) + f(t, YC) At 


Using the sigma notation, the (k + 1)“ value of Y(t) can be written as 
k-i 
Ya+kAth=yot+ ¥ flat+ndt, Y(a+nAt)) At. 


n=0 


This equation may also be written in the manner of a Riemann sum with b =a +-k At: 


b 
¥(b) = yo + © (GY (0) At. 


In the simple case 


y® =f, 


the Euler approximation is just yg plus the Riemann sum, 
b 
Y(b) = yo + DY S(N AL 


which is approximately equal to y9 plus the integral 


b 
yb) = yo + [ fat. 


EXAMPLE 1 Compute the Euler approximation to the initial value problem 
y’ =t— yr y(0) = 0 
forO <t < 1, with Ar = 0.2, 


Notice that the differential equation is not linear because of the y?, and the 
variables are not separable, so we cannot solve the equation by the methods 
of the preceding sections. Given Y(f), the next value Y(t + At) is computed 
by the rule 


Y(t + At) = Y(t) + (t — ¥(t)*) At. 
We record the values in a table. The third column gives the change in Y(t). 


The graph of Y(t), shown in Figure 14.4.2, is obtained by connecting the 
points (t, Y(t)) in the table by straight lines. 


At = 0.2 
t Y@) Y(t + At) — Y(t) = (t — Y(t)”) At 
0.0 0.0 0.0 
0.2 0.0 0.04 
0.4 0.04 0.0797 
0.6 0.1197 0.1171 
0.8 0.2368 0.1488 
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Figure 14.4.2 Example 1 


Work the same problem with At = 0.1. 


At = 0.1 
t Y(t) Y(t + Ad) — Y(t) = (t — Y(t”) At 

0.0 0.0 0.0 
0.1 0.0 0.01 
0.2 0.01 0.02 
0.3 0.03 0.0299 
0.4 0.0599 0.0396 
0.5 0.0995 0.0490 
0.6 0.1486 0.0578 
0.7 0.2063 0.0657 
0.8 0.2721 0.0726 
0.9 0.3447 0.0781 


We now consider Euler approximations with infinitesimal At. These approxi- 
mations cannot be computed directly but are useful in showing that a differential 
equation has a solution. 

The Euler approximation Y(t) depends on both t and the increment size At. 
Now let At be positive infinitesimal. By the Transfer Principle, Y(t + At) is still 
given by the rule 


Y(t + At) — Y(t) = f(t, YO) At. 


Intuitively, the graph of Y(t) as a function of t is formed from infinitesimal line 
segments, and the segment from f to t + At has slope f(t, Y(‘)), as in Figure 14.4.3. 

The next theorem shows that the Euler approximation for infinitesimal At 
is infinitely close to a solution of the initial value problem. 


EXISTENCE THEOREM 


Let At be positive infinitesimal and let Y(t) be the Euler approximation of 
the initial value problem 


(1) Y=fty, YWO=Yo 
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\ 


| 


Figure 14.4.3 


with increment At. Let b be a real number greater than a, and suppose that 
Y(0) is finite for all t between a and b. Then for real numbers t in the interval 
La, b], the function y(t) given by 


y(t) = se(Y() 


is a solution of the initial value problem (1). 


Discussion The theorem shows that the initial value problem (1) has a solution 
as long as Y(t) remains finite. The solution is found by taking the standard 
part of Y(t). When Y(t) becomes infinite, we say that an explosion occurs 
(see Example 2 in Section 14.1 and Example 3 in this section). 


PROOF OF THE EXISTENCE THEOREM At t =a, \(a) = st(Y(a)) = yo. Let M 
be the largest value of | f(t, Y(t))| for 4, a partition point between a and b, 
Then M is finite. Since Y(t) never changes by more than M At from one 
partition point to the next, we always have 


| Ya) -— Y(s)| < M|t — s|. 
Taking standard parts, we see that for real s and ¢ in the interval [a, b], 
Iy(t) — y(s)| s M|t — s|. 


By the Transfer Principle, this also holds for all hyperreal s and t between 
aand b. Then for anya <t <b, 


Y() ® Y(st()) & y(st(t)) * yo) 
and hence, because f(t, z) is continuous in z, 
FG YO) & FG yO). 
Let h(t) be the real function 
A(t) = f(t, x). 


Since Y is an Euler approximation, 
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YQ)=yYot+ DY f(s, Y(s)) At 
for each real point t between a and b. But A(s) = f(s, y(s)) = f(s, Y(s)), so 
Y(Q) ® yo + ¥ h(s) At. 


This is just the Riemann sum of h. Taking standard parts, we get the integral 
of h: 


0) =o 4 { h(s) ds. 


Finally, by the Fundamental Theorem of Calculus, 


Y(t) = hO =fG yO). 


Thus y(t) is a solution of (1) as required. 


To apply the Existence Theorem, we need a way of checking that Y(t) is 
finite. Here is a convenient criterion. 


LEMMA 


Let Y(t) be an Euler approximation of the initial value problem (1) with 
infinitesimal At, and let M and b be finite. 


@ Iftft,y)| <M for alla<t <b and all y, then Y(t) is finite for all 
a<t<b. 

Gi) If | f(t, y)| < M for all a<t <b and all y within M +(t — a) of yo, 
then Y(t) is finite for alla <t < b. 


PROOF (i) Since Y(t) cannot change by more than M At from one partition point 
to the next, we have 


\YO-—yol <M-(t-—a)<M-(b—- a). 
M- (b — a) is finite, so Y(t) is finite. 
The proof of (ii) is similar. 


Discussion The lemma is illustrated in Figure 14.4.4. Choose a positive real number 
M. Part (i) of the lemma says that if f(t, y) is between — M and M, everywhere 
in the vertical strip between t = a and t = b, then Y(t) is finite fora <t <b. 
Part (ii) of the lemma says that if f(£, y) is between —M and M, everywhere 
in the shaded triangle, then Y(t) is finite for a < t < b. Part (ii) is stronger 
because the shaded triangle is a subset of the vertical strip. The proof shows 
that Y(t) stays within the shaded rectangle fora <t < b. 


The lemma and the Existence Theorem combined show that if we can find 
an M such that f(t, y) is continuous and f(t, y) is between —M and M every- 
where in the shaded triangle, then the initial value problem (1) has a solution 
y(t) for a < t < b. The proof also shows that the solution y(t) is within the 
shaded triangle for a < t < b. 
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Figure 14.4.4 


EXAMPLE 2 Show that the initial value problem 
y=t—y%,  yQ)=0 
from Example 1 has a solution for0 <t < 1. 
Let At be infinitesimal and form the Euler approximation Y(¢) with increment 
At. Apply the lemma with M = 1, b = 1. In this example, 
Ifeyl=l?-th<1 
whenever 0 <t <1 and —1< y<1. Therefore, by the lemma, Y(t) is 


finite for 0 < t < 1. By the Existence Theorem, the standard part of Y(t) is a 
solution forO <t <1. 


Here is another theorem that shows that in most cases the solution is unique 
and is close to the Euler approximations for small real increments At. In this theorem, 
we shall write Y,,(f) instead of Y(t) to keep track of the fact that Y(t) depends on At 
as well as on f¢. 


UNIQUENESS THEOREM 


Assume the hypotheses of the Existence Theorem and also that f(t, y) is 
smooth; that is, the partial derivatives of f are continuous. Then the initial 
value problem (1) has only one solution y(t) for t in [a,b]. Furthermore, the 
Euler approximations Y,(t) approach y(t) as the real number At approaches 
zero; that is, 


lim, Ya = xO 


for each t in [a, 6]. 


14.4 EXISTENCE AND APPROXIMATION OF SOLUTIONS 


We shall not give the proof. The Uniqueness Theorem tells us two important 
things about differential equations in which f(t, y) is smooth. First, it tells us that a 
particular solution of such a differential equation will depend only on the initial 
condition. Thus if an experiment is accurately described by a differential equation 
with f(t, y) smooth, then repeated trials of the experiment with the same initial 
condition will give the same outcome. Second, it tells us that the Euler approximations 
will approach the solution of the differential equation as At approaches zero. Thus 
we can get better and better approximations of the solution by taking Ar small. 


EXAMPLE 2 (Continued) The function f(t, y) = t — y? is smooth. The Uniqueness 
Theorem shows that the initial value problem of Example 1 has just one 
solution y(t) for 0 <t <1. Moreover, the Euler approximations Y(t) get 
close to y(t) as the real increment At approaches zero. Thus the approxima- 
tions computed in Example 1 really are approaching the solution. 


We conclude with an example of an explosion and an example with more 
than one solution. 
EXAMPLE 3 (An Explosion) The initial value problem 
y=y,  y0=1 


may be solved by separation of variables: 


y? dy = dt, 
= 
—. =t+C, 
y(t) 
wt) = ~(t+ 0)" (general solution), 
i2-O2es). Cast 
yt) = -—17! (particular solution). 


The graph, shown in Figure 14.4.5, approaches infinity as t approaches 1 
from the left. The function y(t) = (1 — t)~' is a solution for 0<t <1, 
and the solution has an explosion at t = 1. 


y 


~ 


Figure 14.4.5 Example 3 


871 


872 14 DIFFERENTIAL EQUATIONS 


The Euler approximation Y(t) with a real increment Ar can be computed 
even for ¢ greater than 1, but will approach y(t) only for 0 <t < 1. For 
infinitesimal Az, the Euler approximation will be finite and infinitely close 
to y(t) when ¢ is in the real interval [0, 1). Y(2) will keep on increasing and 
will be infinite for all ¢ with standard part > 1. 


Continuing the example, compute the Euler approximation Y(t) for At = 0.2 
and 0 < ¢ < 2, and compare the values with the solution y(t) = (1 — N73 
for 0 <t <1. The results are shown in the next table and are graphed in 
Figure 14.4.5. 


At = 0.2 
t y(t) Y(t) Y(t)? At 
0.0 1.0 1.0 0.2 
0.2 1.25 1.2 0.288 
0.4 1.6667 1.4488 0.4428 
0.6 2.5 1.9309 0.7456 
0.8 5.0 2.6764 1.4327 
1.0 ao 4.1091 3.3770 
1.2 7.4861 11.2084 
1.4 18.6945 69.8966 
1.6 88.5910 1569.6736 
1.8 1658.26 549968.3 
2.0 551627.6 


EXAMPLE 4 (Nonuniqueness) The initial value problem 
y =3y, (0) = 0 


has infinitely many solutions. The graphs split apart, as shown in Figure 
14.4.6. 


Figure 14.4.6 Example 4 
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One solution is the constant solution y(t) = 0. A second solution is found 
by separating variables: 
ay 7° dy = dt, 
ylBuat+c, C=0 
y=. 
This solution can be checked by differentiation: 
yi = 3t? = 3(09)79 = 3y?/9. 

The other solutions go along the line y(t) = 0 in some interval and branch 
off the line y(t) = 0 to the right and left of the interval. The full list of solu- 
tions Is: 

(t-—a) fort<a 

y(t) = 40 fora<t<b 

(t—b)> forb<t 
Here a is either a real number or — ©, b is either a real number or + 0, 
and a < 0 < b. In the case that a = —co and b = +0, the solution is the 
constant function y(t) = 0. 
These solutions all have the same initial value y(0) = 0. The Uniqueness 


Theorem does not apply in this example because the function f(t, y) = 3y?? 
has no derivative at y = 0, so that f(£, y) is not smooth at y = 0. 


PROBLEMS FOR SECTION 14.4 


In Problems 1-4, compute the Euler approximation to the given initial value problem with 
At=O.1forO0<t<l1. 


1 
3 
5 


y=tly, y0)=1 2 y=t+y, yO)=1 
y’ = cos(t + y), yO) = 1 4 y =tlny, y(O) = 2 
Show that the initial value problem 

y=14+t’, y0)=0 
has a unique solution y(t) for 0 <t < 4. 
Show that the initial value problem 


F y 
=ttot 
is 2 


has a unique solution y(t) for 0 < t < 4. 
Show that the initial value problem 

y = arctan (ft + e”), y(0) = 4 
has a unique solution y(t) for 0 < t < oo. 
Show that the initial value problem 

y=t+exp(-y’), yO)=1 
has a unique solution y(t) for 0 < t < @. 
Show that the initial value problem 

y=y', (0) =0 

has infinitely many solutions for 0 < tf < oo. 
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10 Show that the initial value problem 
y=e«lt—y)'?, yO)=1 
has infinitely many solutions for0 <t < a. 
11 Suppose that f(t, y) is continuous for all t and y. Prove that for each point (a, yo), the 
initial value problem 
y =cos(f(,y)), Wa) = Yo 
has a solution y(t),a <t < 0. 
12 Suppose that f(t, y) and g(t) are continuous for all f and y and that | f(4 y)| < g(O 
for all t and y. Prove that for each point (a, yo) the initial value problem 
y=fGy), Wa) =o 
has a solution y(t) fora <t < @. 
13 Suppose that f(t, y) is continuous for all ¢ and y. Prove that for each point (a, 9) there 
is a number b > a such that the initial value problem 
Y=fy),  y@) = yo 
has a solution y(t)a <t < b. 
COMPLEX NUMBERS 


This section begins with a review of the complex numbers. Complex numbers are 
useful in the solution of second order differential equations. The starting point is 
the imaginary number i, which is the square root of —1. The complex number system 
is an extension of the real number system that is formed by adding the number i 
and keeping the usual rules for sums and products. The set of complex numbers, or 
complex plane, is the set of all numbers of the form 


z=x4+iy 


where x and y are real numbers. The number x is called the real part of z, and y is 
called the imaginary part of z. A complex number z can be represented by a point 
in the plane, with the real part drawn on the horizontal axis and the imaginary part 
on the vertical axis, as in Figure 14.5.1. The sum of two complex numbers is computed 
in the same way as the sum of two vectors, 


(a+ ib) +(e + id) =(a+c)+ ib + d). 


iy 


Figure 14.5.1 


14.6 COMPLEX NUMBERS 
The product of two complex numbers is computed using the basic rule i? = —1 
and the rules of algebra: 


(a + ib) *(c + id) = ac + ibe + iad + i2bd = (ac — bd) + i(bc + ad). 


EXAMPLE 1 Compute the product of 3 + i6 and 7 —i. 
(3 + i16)-(7 —i) = (3 +7 —6°(—1)) + (6-74 3+(—1)) = 27 4 139. 


The complex conjugate Z of z is formed by changing the sign of the imaginary part of z: 
a+ib=a-— ib. 

The product of a complex number with its conjugate is always a nonnegative real 

number, computed as follows. 

(1) (a + ib)(a — ib) = a? — iab + iab + b? = a? + b?. 


The quotient of two complex numbers can be computed by multiplying the numerator 
and denominator by the conjugate of the denominator, as follows. 


' atib_ (a+ibic—id) _ac+bd  .—ad+be 
| e+id (c+idlec-id) C+@'' O4+d 


EXAMPLE 2 Compute the quotient (1 + i)/(1 — 72). 


1+i (@+i-(+i2) 9 -1—1-2)4+ 10-24 161) 
i=? CD04 D 1 
ioe 1 3 


Sree to 


In the real number system, a positive real number b has two square roots, 
Jb and ~./b, and negative real numbers have no square roots. In the complex 
number system, a negative real number —b has two imaginary square roots, i,/b 
and —i,/b. The quadratic formula gives the roots of any second degree polynomial 
in the complex number system. 


QUADRATIC FORMULA 


The roots of the polynomial 
: az?>+bz+c  wherea#0 
in the complex number system are given by 


_ b+ /b? — 4ac 
2a : 


The number b? — 4ac is called the discriminant. If a, b, and c are real, there 
are three cases: 


Case 7 If b? — 4ac > 0, there are two real roots. 
Case 2 If b? — 4ac = 0, there is one real root. 
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Case 3 If b? — 4ac < 0, there are two complex roots, which are complex conjugates 
of each other. 


EXAMPLE 3 Find the roots of the polynomial z* + z + 2 in the complex number 
system. 


Sb alee: He a7 1 v7 
as 5 = 


It is often useful to represent a complex number in polar form. A point 
(x, y) in the plane has polar coordinates (r,@) where x = rcos 6, y =rsin@. The 
complex number x + iy may be written in the polar form 


x + iy = r(cos 6 + isin @). 


The coordinates r and 8 can always be chosen so that r > 0 and —z < @ <7. The 
number r, which is the distance of the point (x, y) from the origin, is the absolute value 
of the complex number x + iy: 


r= [x + iy| = (x? + y?)"?. 


The formula (1) for the product of a complex number and its conjugate may now 
be written in the short form 


zz= |z2/. 
The real number @ is an angle in radians and is called the argument of x + iy. The 
argument can be computed by using the formula 


Yi 
tan~ = 60, 
x 


and then choosing @ in the correct quadrant. The polar form of a complex number 
is illustrated in Figure 14.5.2. This figure is sometimes called the Argand diagram 
of the complex number. The complex number with absolute value one and argument 0 
is sometimes called cis 0: 


cis@ = cos @ + isin @. 


Figure 14.5.2 
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In Figure 14.5.3, we see that in the complex plane, cis @ is on the unit circle at an 
angle @ counterclockwise from the x-axis. Using the symbol cis @, the polar form 
can be written 

x+iy=rcis 0. 


cis 0 


Figure 14.5.3 


EXAMPLE 4 Write the complex number z = —2 + i2 in polar form. 


The absolute value of z is |z] = (2? + (—2)?)!? = JB. To find the argument 
9, we use tan 6 = 2/(—2) = —1. Since zis in the second quadrant (x negative 
and y positive), 86 must be 32/4. Thus 


z= /Bcis 


cis @ is helpful in computing products, quotients, and powers of complex 
numbers. Using the addition formulas for sines and cosines, we can prove the product 
formula 


(2) (r cis 8) + (s cis @) = rs cis (0 + ¢). 


In words, this formula states: To multiply two complex numbers, multiply the absolute 
values and add the arguments. There is a similar formula for quotients: 
rcis@ or 


(3) = -—cis (0 — @). 


scish 8 


To divide two complex numbers, divide the absolute values and subtract the arguments. 


EXAMPLE 5 Using the polar form, find the quotient (1 + )/( — ?). 


In polar form, 


L+i=V2eiss, 1 i= /es(-3), 
1+i SS 

ti V2 o(m _ : = cis(= 

1—-i J2 4 4 2 
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Using the product formula (2) n times, we get a formula for the n™ power of a complex 
number, 
(4) ( cis 8)" = r" cis (n€). 
This formula in the case r = 1 is called De Moivre’s Formula, 
(cos @ + isin 6)" = cos (n8) + isin (18). 
We can see from the power formula (4) that the complex number r cis @ has the 
square root Jr cis (8/2). In fact, each complex number except zero has two square 
roots, 
7 1/2 : 8 
(5) (r cis #)* = +./r cis 5. 
EXAMPLE 6 Find the square roots of i. 
By the computation in Example 3, the polar form of i is i = cis (1/2). 
pe mp Teis = + (cos + isin) 


(2.02) 


The two square roots of i are shown in Figure 14.5.4. 


VE(1 + i) 


2 +8 


Figure 14.5.4 


We now turn to complex exponents, which are useful in the study of differ- 
ential equations. In order to give a meaning to an exponent e’, we consider infinite 
series of complex numbers. The sum of an infinite series of complex numbers is 
defined by summing the real and imaginary parts separately. If z, = x, + iy,, and 
the series }’ x, and }° y, both converge, the sum of the series }" z, is defined by the 
formula 


ao 


XY Z, = oy X, +i ms Vue 
n=0 n=0 


n=0 


In Chapter 9, we found that for real numbers z the exponent e’ is given by the power 
series 
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z? a zt 


When z is a complex number, this formula is taken as the definition of e’. It can be 
shown that the power series converges for every z and that the exponential rule 
e+? = e%e? holds for complex exponents. In the case that z is a purely imaginary 
number z = iy, the power series takes the form 


Using the power series for cos y and sin y, we obtain Euler’s Formula: 
e” =cosy + isin y = cis y. 
When z is a complex number z = x + iy, the exponent e’ is given by the formula 


ext? = eel” = e*(cos y + isin y). 


EXAMPLE 7 Find e~ 2+ '#/3. 


fas 
: e 2 tials — | cos + isin] = (5 + 2) , 


) 


In Chapter 8, the hyperbolic cosine and hyperbolic sine were defined in 


terms of e* by the equations 
e+te* e—e* 
cosh x = —_~—_,, sinh x = 
2 2 


Euler’s Formula leads to similar equations for the cosine and sine. 


e” = cos y + isin y, e ” = cosy — isin y, 
Jy: 
e% +e"? : ev — eo 
cos y = ———— sin y = ——_-. 
a 4 a i2 


In the next section we will make use of complex valued functions, that is, 
functions f(t) that assign a complex number z = f(t) to each real number t. The 
derivative of a complex valued function is obtained by differentiating the real and 
complex parts separately. Thus, if h(t) = f(t) + ig(t), where g and # are real functions, 
then h'(t) = f’@® + ig’. 

For example, if h(t) = e”, where r = a + ib is a complex constant, then 

h(t) = e“(cos bt + isin bt), 
h'(t) = ae“(cos bt + isin bt) + be“(—sin bt + icos bt) 
= (a + ib)e“(cos bt + isin bt) =re™. 

Summing up, the usual rule (e")’ = re” still holds when r is a complex 
constant. We can also consider complex valued differential equations. The example 
we shall need is the homogeneous linear differential equation 


z+rz=0, 
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where r is a complex constant. The general solution of this equation is 
2(t) = Ce", 


where C is a complex constant. This solution can be checked by differentiation as 
before. 


EXAMPLE 8 Solve the complex initial value problem 
z+(34+ i4)z=0, 2(0) = et F!?, 
The general solution is 
2(t) = Ce~ 3+ i4H, 


Substituting the initial value at t = 0, e!*!* = C. The particular solution 
is then 


2(t) = e! Fi26- GB tidy = e} —3i+ (2 — 41) 
The solution may also be written in polar form using Euler’s Formula, 


2(t) = e! 3 cis (2 — 41). 


PROBLEMS FOR SECTION 14.5 


In Problems 1-6, put the complex number in the form a + ib. 


1 (7 — 14) + G + i2) 2 (4 — i6) — (8 — i) 
3 (—4 + 12) -(1 — i2) 4 (3 + i) *(—2 — 16) 
5 (1 — i238 +d 6 (7 + i3)/Q2 — i5) 
In Problems 7-12, find the roots of the given equation. 

7 2? —8z+ 16=0 8 z4+6z+9=0 
9 2+4+25=0 10 z7+ 100=0 

11 z+2z24+5=0 12 2+z+3=0 


In Problems 13-18, put the complex number into the polar form r cis 0. 
13 iS 14 —i3 
15 —3 — 73 16 —4+4+ i4 


17 esi 18 24 12/3 


In Problems 19-24, use the polar form to simplify the given expression. 


19 (2 + i2)(—3 + 13) 20 (—4 — i4)(5 — iS) 

21 (/3 + yd —) 22 (+ a/ — 1/3) 

23 (i — i) 24 (1/3) 

In Problems 25-28, compute both square roots of the given complex number using the polar form. 
25 1+i 26 —1+i 

27 —i4 28 agen i) 


In Problems 29-32, put the given exponent in the form a + ib. 


29 ew 3 tin 30 e2—in 
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el inl4 32 e ini6 


In Problems 33-36, find the derivative. 


33 
35 


z= ef By 34 z= elt ti2y 


z= ef2 IN +G + i2)t 36 z= er tiitian 


In Problems 37-40, find the general solution of the complex differential equation. 


37 
39 


z+ (2 —-i3)z=0 38 z — i4z=0 
z+ (—3 4+ i5)z =0 40 2 +(-2-i)z=0 


In Problems 41-44, solve the given complex initial value problem. 


41 
42 
43 
44 
45 
46 
47 


48 


49 


50 


51 


52 


53 


zZ+(-24+i)z=0, 20)=e! 

zZ+3+i4z=0, 2(0)=e*** 

z+(-2+iz=0, z(0) = 4 

z+(3+4 i4z =0, z(0) = —1 

Show that for any complex number z, z + Z is a real number. 
Prove that the conjugate of a complex number r cis 6 is r cis (— 8). 


Prove that for any nonzero complex number z, 2/Z = cis (20), where 6 is the argument 
of z. 


Prove that for any two complex numbers u and z, the sum of the conjugates of u and z 
is equal to the conjugate of the sum of u and z, and similarly for products. In symbols, 


@+Z=ut+Z  (i)+@=uez. 
Prove that for any two complex numbers u and z, 
lu + z| < |ul + |Z]. 
Use. De Moivre’s Formula with n = 2, 
cos (20) + isin (20) = (cos 6 + isin 6)?, 
to obtain expressions for cos (20) and sin (26) in terms of cos @ and sin 6. 
Use De Moivre’s Formula with n = 3, 
cos (36) + isin (30) = (cos @ + isin 0)°, 
to obtain expressions for cos (30) and sin (30) in terms of cos @ and sin 0. 
Find the solution of the initial value problem 
z+ (a+ ib)z =0, 2(0) = e°*i4 
where a, b, c, and d are real numbers. 


Show that every solution of the differential equation z’ + ibz = 0 has constant absolute 
value (where b is a real number). 


SECOND ORDER HOMOGENEOUS LINEAR EQUATIONS 


A second order differential equation is an equation involving an independent variable t, 
a dependent variable y, and the first two derivatives y’ and y”. The general solution 
of a second order differential equation will usually involve two constants. Two initial 
values are needed to determine a particular solution. A second order initial value 
problem is a second order differential equation together with initial values for y and y’. 
This section gives a solution method for second order equations of the following 
very simple type. 
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SECOND ORDER HOMOGENEOUS LINEAR DIFFERENTIAL EQUATION 

WITH CONSTANT COEFFICIENTS 

(1) ay” + by’ + cy =0, 
where a, b, and c are real constants and a,c # 0. 

Discussion If a = 0, the equation is a first order differential equation by’ + cy = 0. 
If c = 0, the change of variables u = y’ turns the given equation (1) into a 
first order differential equation au’ + bu = 0. In each of these cases, the 


equation can be solved by the method of Section 14.1 or 14.2. (After finding 
u, y can be found by integration because y’ = u.) 


In Section 14.2 we found that the first order homogeneous linear differential 
equation with constant coefficients, y’ + cy = 0, has the solution y = e~“. To get 
an idea of what to expect in the second order case, Jet us try to find a solution of 
equation (1) of the form » = e" where r is a constant. Differentiating and substituting 
into equation (1), we see that 


a(e")" + b(e"y + ce = arre" + bre” + ce" 
= (ar* + br + ce". 
This shows that y(t) = e” is a solution of equation (1) if and only if 
ar? +br+c=0. 


We should therefore expect that the solutions of the equation (1) will be built up 
from the functions y(t) = e” where r is a root of the polynomial az* + bz + c. We 
shall state the rule for finding the general solution of the equation (1) now and 
prove the rule later on. 


METHOD FOR SOLVING A SECOND ORDER HOMOGENEOUS DIFFERENTIAL 
EQUATION WITH CONSTANT COEFFICIENTS 


(1) ay” + by’ +cy =), a # 0. 


Step 7. Form the characteristic polynomial 


az? + bz+. 
Find its roots by using the quadratic equation or by factoring. 

Step 2. The general solution is described by three cases. 
Case 7 Two distinct real roots: 2.7 2 = "s; 

y = Ae™ + Be™. 
Case 2 One real root: z=, 

y = Ae™ + Bte”™. 
Case 3 Two complex conjugate roots: z=art ip. 

y = e“[A cos (Bt) + Bsin (f2)]. 


Step 3 If initial values for y and y' are given, solve for A and B, and substitute to 
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obtain the particular solution. The two initial values will specify the position 
and velocity at one time: 


Y= Vo and y =o att = fo. 
Discussion The general solution in Case 3 is sometimes written in the same form 
as Case 1 by using complex exponents, 
y = Ce" + De", 
where r=a+ ip, s=o — if, 
C = 4(A — iB), D=4(A + iB). 
To show that the two forms of the solution are really the same, use the 
complex exponent formula 


e** 8 — e%(cos B + isin fp) 


from the preceding section. 


EXAMPLE 1 Find the general solution of 
y"—-ow’y=0, wF0. 


Step 1 The characteristic polynomial is z? — w. It has two real roots, z = and 
Z= —O. 


Step 2 The general solution is 
y = Ae™ + Be™™, 


where @ is constant. 


EXAMPLE 2 Find the solution of the initial value problem 
y’-y—2y=0, yO)=5, y(0)=0. 


Step 1 The characteristic polynomial z* — z — 2 has two real roots, z= —1 
and z = 2. 


Step 2. The general solution is 
y = Ae™' + Be*, 
Step 3 The initial value y(0) = 5 gives the equation 
S=A+B. 


To get a second equation, we differentiate the general solution and substitute 
the initial value for y’(0). 


y’ = —Ae™' + 2Be”', 0= —A + 2B. 


The solution of the two equations for A and B is 
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The particular solution of the initial value problem, shown in Figure 14.6.1, is 
10 5 
y= (Fe + (5)e* 


(0, 5) 


(1, 0) 


Figure 14.6.1 Example 2 


EXAMPLE 3. Find the solution of the initial value problem 
y"— 4y + 49 =0, (0) = —3, y'(0) = 1. 
Step 1 The characteristic polynomial z* — 4z + 4 has one real root, z = 2. 
Step 2 The general solution is 
y = Ae?! + Bre, 
Step 3 Substitute 0 for t and —3 for y. 
—3 = Ae® + B+0 +e, A= —3, 
Compute y’ for the general solution. 
y’ = 2Ae”! + 2Bte* + Be", 
Substitute 0 for t and 1 for y’. 
1 = 2Ae° + 2B *0+e° + Be? = 2A 4B, B=7, 
The particular solution, shown in Figure 14.6.2, is 


y = —3e7! + Tte*. 


EXAMPLE 4 Find the solution of 
2y" + 18y = 0, y(0) = 2, y'(0) = 15. 
Step 1 The characteristic polynomial is 2z? + 18, and its roots are z = +13. 
Step 2. The general solution is 
y = Acos (3t) + Bsin (30). 
Step 3 Substitute 0 for t and 2 for y. 
2=Acos0+ Bsin0O = A. 


14.6 SECOND ORDER HOMOGENEOUS LINEAR EQUATIONS 885 


~ 


(0, —3) 


Figure 14.6.2 Example 3 


Compute y’(t) for the general solution. 
y = —3A sin Gt) + 3Bcos (3t). 
Substitute 0 for t and 15 for y’. 
15 = —3Asin0 + 3BcosO = —3+0+3°B, B=5. 
The particular solution, shown in Figure 14.6.3, is 
y = 2cos (3t) + 5 sin (32). 


(0,2) 


Figure 14.6.3 Example 4 


GRAPHS OF SOLUTIONS OF SECOND ORDER HOMOGENEOUS LINEAR EQUATIONS 


Our next topic is graphs of solutions of second order homogeneous linear equations. 
Several cases arise, including simple, damped, and overdamped oscillations. 
Consider the second order homogeneous linear differential equation 


(1) ay” + by’ +c=0. 
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We shall concentrate on the case that 
a>QO, b>0, c>0. 


(We can always make a positive by changing all the signs if a is negative. The cases 
with negative b or c are considered in the problem set at the end of the section.) 
The equation has the characteristic polynomial 


az? + bz+c. 
Let d be the discriminant of this polynomial, 
d= b? — 4ac. 


Simple Oscillation This type of solution arises when b = 0. Since a and c are assumed 
to be positive, the discriminant d = —4ac is negative, and the characteristic poly- 
nomial has two purely imaginary roots, +if. The general solution of equation (1) is 


y(t) = Acos (pt) + Bsin (ft). 


It is helpful to put this equation in a different form. The point (A, B) is on the circle 
with center at the origin and radius C = (A? + B?)'/?, There is thus an angle @ 
for which 


A =Ccos 8, B=Csin8, 
as in Figure 14.6.4. The angle @ can be computed as follows: 
B Csin@ — 


A Ccos0 


B 
tan 0, 8 = arctan () ; 


(A, B) 


B=Csin@ 


Figure 14.6.4 


Using the formula for the cosine of the difference of two angles, cos(¢@ — @) = 
cos (#) cos (8) + sin (f) sin (@), we find that 


y(t) = C cos (Bt) cos (0) + C sin (Br) sin (@) = C cos (Bt — 8), 
so that 
y(t) = C cos (ft — 8). 


The number C is called the amplitude, because the cosine curve oscillates between 
C and —C. The number f is called the frequency, because the curve will complete 
B cycles each 2x units of time. The number 27/f is called the period, because each 
cycle is 2z/f units long. The angle 6 is called the phase shift. Thus the graph of each 
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particular solution is a cosine wave with amplitude C, period 27/f, and phase shift 6, 
as illustrated in Figure 14.6.5. 


y 


™ 


Figure 14.6.5 


Damped Oscillation This case arises when b is positive and the discriminant is 
negative, b > 0 and d < 0. Here b is in the range 0 < b < ,/4ac. The roots of the 
characteristic polynomial are complex conjugates, « + if. The real part « = —b/2a 
is negative. The general solution is 


y(t) = e“[A cos (Bt) + B sin (Bt)]. 


Each particular solution oscillates with period 2x/f, but the amplitude dies down 
exponentially as t increases, as in Figure 14.6.6. 


Figure 14.6.6 


As in the case of simple oscillation, the solution may be written in the form 
y(t) = e“C cos (Bt — 8), 


where C = (A? + B?)!/? and @ is a constant angle. The amplitude at time ¢ will then 
be e“'C, which is decreasing because « is negative. 


Critical Damping This case arises when b is positive and the discriminant is zero, 
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so’ that b= /4ac. The characteristic polynomial has one negative real root, 
r = —b/2a. The general solution is 

y(t) = Ae" + Bte™. 


Each particular solution will approach 0 as t approaches infinity and will never 
complete one oscillation. A solution can cross the x-axis once, but never more than 
once. See Figure 14.6.7. 


y 


e 7 


Figure 14.6.7 


Overdamping This case arises when 5 is positive and the discriminant is positive, 
so that b > ./4ac. The characteristic polynomial has two real roots, r and s. Since 
a, b, and ¢ are positive, the characteristic polynomial cannot have any positive or 
zero roots. Therefore both roots r and s are negative. The general solution is 


y(t) = Ae™ + Be™. 


Again, each particular solution approaches zero as t approaches infinity and will 
never complete one oscillation, as in Figure 14.6.7. 

The differential equations of this section provide simple models for a variety 
of physical systems that oscillate, such as mass-spring systems and electrical networks. 
When a horizontal spring of natural length L is compressed a distance x, it exerts a 
force of approximately F = —kx. k is called the spring constant and depends on the 
particular spring. The negative sign indicates that the force is in the opposite direction 
from x, as in Figure 14.6.8. A mass m is attached to the end of the spring. From 
Newton’s Law, 

F = ma = mx", 


we obtain a second order differential equation for the position x(t), 
mx" = —kx, 


or mx" + kx = 0. 


Figure 14.6.8 
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Both constants k and m are positive, and the solution is the simple oscillation 
x(t) = A cos (ft) + B sin (ft) 


where B = (k/m)"?. 

A mass-spring system immersed in oil (such as an automobile shock 
absorber) is subject to a damping force bx’(t), which is proportional to the velocity 
x(t) but in the opposite direction. This additional force will slow down the motion 
of the spring and lead to a damped oscillation. The force is approximately 


F= —bx’ — kx 
and thus satisfies the differential equation 
mx” + bx’ + kx = 0. 


When the damping constant b is between 0 and (4mk)!/”, the solution will be a damped 
oscillation. The greater the value of b, the more quickly the oscillation will be damped 
down. When b is equal to (4mk)"/?, the solution will be critically damped; when b 
is greater than (4mk)‘/?, the solution will be overdamped. 


EXAMPLE 5 Suppose a mass-spring system 
mx” + bx’ + kx 


has spring constant k = 5, damping constant b = 4, and mass m = 1. At 
time t = 0, the position is x(0) = 1 and the velocity is x‘(0) = 2. Find the 
position x(t) as a function of time. 


The differential equation is 
x" + 4x' + 5x = 0. 


Step 1 The characteristic polynomial z* + 4z + 5 has roots —2 —iand —2 +i. 
(These roots can be found using the quadratic equation.) 


Step 2. The general solution is 
x(t) = e*[Acost + Bsint]. 
Step 3 Find A and B using the given initial values. 
1= e°[Acos0 + Bsin 0], A=l1. 
Compute x’(t) and substitute to find B. 
x(t) = —2e7*[A cost + Bsint] + e *[—Asint + Bcost], 
2 = —2e°[cos 0 + Bsin 0] + e°[—sin0 + Bcos 0], 
B=4., 
The particular solution is 
x(t) = e” [cost + 4sin ¢]. 


This is a damped oscillation with period 2z. 


Let us now justify the solution method given at the beginning of this section. 
We may take the coefficient of y” to be one and consider the second order differential 
equation 
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(2) y’ + by’ +cy =0. 
Let r and s be the roots of the characteristic polynomial, so that 
2+obzt+eos(z—r(z—s)= 2? —1z—sz+ 4s. 


If + and s are distinct (either two different real numbers or complex conjugates), we 
must show that the general solution is 


y = Ae" + Be™. 
If » = s, we must show that the general solution is 
y = Ae™ + Bte™. 


The plan is to break the second order equation (2) into a pair of first order differential 
equations. It is useful to use the symbol D for the first derivative and D? for the second 
derivative with respect to t. Thus 


Dy = y’, Dy = yp”, 
The differential equation (2) can then be written in the form 
(3) (D? + bD + c)y = 0. 


We now wish to “factor” the expression D* + bD +c as if it were a polynomial. 
It is to be understood that (D — r)(D — s)y means (D — r)u where uw is the function 
(D — s)y = y’' — sy. Thus, using the Sum and Product Rules for derivatives, 


(D — r)(D — s)y = (y’ — sy)’ — rQ’ — sy) = y" — sy’ — ry’ + rsy 
= y” + by’ + cy = (D? + bD + oy. 


This shows that the second order equation (3) is equivalent to the pair of first order 
equations 


(4) (D — ru 
(5) (D —s)y 
Equation (4) is a homogeneous linear equation whose general solution is 


u= Ke", 


0, 


ur 


Equation (5) may now be put tn the form 
(6) y —sy = Ke", 
This first order linear equation was solved in Example 3 in Section 14.3. In the case 


r % s, the general solution came out to be 


e+ Be™, 


i= 


Putting A = K/(r — s), we get the general solution 
y = Ae + Be", 
In case r = s, the general solution is 
y = Ate“ + Be*, 


where A = K. In each case we have the required formula for the general solution 
of the original equation (2). 
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PROBLEMS FOR SECTION 14.6 


In Problems 1-8, find the general solution of the given differential equation. 


1 


3 
5 
7 


y+ y —by=0 2 y’ ~4y + 3y=0 
y’ — 10y’ + 25y =0 4 y" +2/2y' + 2y =0 
y" + 16y =0 6 y+ 2y=0 

y” — 2y'+2y =0 8 y”" —6y+13=0 


In Problems 9-16, find the particular solution of the initial value problem. 


9 
10 
11 
12 
13 
14 
15 
16 


yt 6y +5y=0, ypO)=1, yO)=0 
y’—y—ly=0, yO)=0, y(0)=14 

y’ + 12y'4+36y=0, y0)=5, y(0)= —10 
y’—8y +16y=0, »0)=-3, y(0)=4 
yt+5y=0, yO)=-2, yO)=5 
y+y=0,  yn/sy=0,  y'(a/4) = 2 

y+ 12y'+37y=0, yO=4 y'(O=0 
y’+6y + 18y=0, yO)=0, yO=6 


In Problems 17-20, solve the initial value problem and find the amplitude, frequency, and phase 
shift of the solution. 


17 
18 
19 
20 
21 


22 


23 


O 24 


O 25 


yt4y=0, yO=/3, yO=2 

y” + 100y = 0, y(0) = 5, yO) = 50 

y+2y+1l0y=0, yO)=1, yO=1 

y’—8y + 25y=0, 30)=3, yO)=0 

A mass-spring system mx” + bx’ + kx =0 has spring constant k = 29, damping 
constant 6 = 4, and mass m = 1. At time ¢ = 0, the position is x(0) = 2 and the velocity 
is x’(0) = 1. Find the position x(t) as a function of time. 

A mass-spring system mx” + bx’ + kx =0 has spring constant k = 24, damping 
constant b = 12, and mass m = 3. At time t = 0, the position is x(0) = Oand the velocity 
is x’(0) = —1. Find the position x(z) as a function of time. 

Show that if y(t) is a solution of a differential equation ay” + by’ + cy = 0, such that 
y(t) = Oand y(t) = 0 at some time fo, then y(t) = 0 for all t. 

In the differential equation ay” + by’ + cy = 0, suppose that a is positive and c is 
negative. Show that the characteristic equation has one positive real root and one 
negative real root, so that the general solution has the form y = Ae’ + Be where r 
is positive and s is negative. 

In the differential equation ay” + by’ + cy = 0, suppose that a and ¢ are positive 
and b is negative. Show that there are three cases for the general solution, depending 
on the sign of the discriminant d: 


Case 7 If d is positive, the general solution has the form y = Ae™ + Be™ where r 
and s are positive. 

Case 2 If d is zero, the general solution has the form y = Ae” + Bte™ where r is 
positive. 

Case 3 Ifdisnegative, the general solution has the form y = e“'[A cos (Bt) + B sin (Bt)] 
where « is positive, so that the graph is an oscillation whose amplitude is 
increasing instead of decreasing. 
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14.7 SECOND ORDER LINEAR EQUATIONS 


This section contains a method for solving nonhomogeneous second order differential 
equations. As in the previous section, we deal only with linear equations with constant 
coefficients. We consider equations of the following type. 


SECOND ORDER LINEAR DIFFERENTIAL EQUATION WITH CONSTANT 
COEFFICIENTS 


(1) ay” + by’ + cy = f(t) 


where a, b, and c are real constants and a,c # 0. 


A differential equation of this form describes a mass-spring system where an outside 
force f(t) is applied to the mass. The function f(r) is called the forcing term. 

As before, if a = 0 the equation is a first order linear differential equation 
in y, and if c = 0 it is a first order linear differential equation in y’. In each of these 
cases, the equation should be solved by the methods of Section 14.3. Hereafter we 
assume a,c # 0. 

To get started, let us review the first theorem on first order linear differential 
equations. Theorem | in Section 14.3 states that the general solution of a first order 
linear differential equation is the sum 


y(t) + Bx(t), 


where y(t) is a particular solution of the given equation and x(t) is a particular . 
solution of the corresponding homogeneous equation. Here is a similar theorem 
for second order equations. 


THEOREM 1 
Suppose that y(t) is a particular solution of the second order linear differential 
equation 

(1) ay" + by’ + cy =f(), 


and Ax,(t) + Bx,(t) is the general solution of the corresponding homo- 
geneous linear differential equation 


(2) ax” + bx' + cx = 0. 
Then the general solution of the original equation (1) is 
y(t) + Ax, (t) + Bx(t). 


As in the first order case, this theorem is proved using the Principle of 
Superposition. 


PRINCIPLE OF SUPERPOSITION (Second Order) 


Suppose x(t) and y(t) are solutions of the two second order linear differential 
equations 


ax” + bx’ + cx = f(t), 
ay” + by’ + cy = g(t). 
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Then for any constants A and B, the function 
u= Ax + By 
is a solution of the linear differential equation 


au” + bu' +-cu = Af(t) + Bg(t). 


Theorem 1 breaks the problem of finding the general solution of the equation 
(1) into two simpler problems. 


First problem: Find the general solution of the corresponding linear homogeneous 
equation 


ax” + bx'+c=0. 
Second problem: Find some particular solution of the given equation 


ay" + by’ + cy =f (0). 


The first problem was solved in the preceding section. We now present a method for 
solving the second problem. This method is sometimes called the method of judicious 
guessing, or the method of undetermined coefficients. The method works only when 
the forcing term f(t) is a fairly simple function, of the form 


(3) p(t)e™ cos (Bt) + q(te™ sin (fr), 


where p(t) and q(t) are polynomials. However, when it works it is a very efficient 
method of solution. Often f(t) will be of an even simpler form, such as a polynomial 
alone, or a single exponential or trigonometric function. In the case of a homogeneous 
equation, where f(t) = 0, the zero function y(t) = 0 is a particular solution. The 
idea for solving a linear equation is to guess that the differential equation has a 
particular solution, which looks like the forcing term f(t) but has different constant 
coefficients. By working backwards, it is possible to find the unknown constants 
and discover a particular solution. We illustrate the method with several examples. 


EXAMPLE 1 A mass of one gram is suspended from a vertical spring with spring 
constant k = 100, as in Figure 14.7.1. At time t = 0, the mass is at position 
y(O) = 2.cm and has velocity y’(0) = 50 cm/sec. Find the equation of motion 
of the mass. It is understood that there is no damping, and the origin is at 
the point where the spring is at its natural length. 


Figure 14.7.1 


In this problem there are two forces, the force of the spring and the force of 
gravity. The force of gravity is a constant and is equal to mg dynes, always 
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in the downward direction. The system is described by a second order linear 
differential equation with a constant forcing term, 


my” + ky = —mg. 
In this case k = 100, m = 1, and g = 980, so the differential equation is 
y” + 100y = — 980, 


To solve the problem, we first find some particular solution of the differential 
equation (4), then use Theorem 1 to find the general solution, and finally 
substitute to find the particular solution for the given initial values y(0) = 2 
and y’(0) = 50. 


Since the forcing term is a constant f(t) = —980, we guess that the differ- 
ential equation (4) has a particular solution, which is a constant. In this 
example it is easy to see by inspection that the constant function 
980 
th) = —— = —98 
WO) = 95 


is a particular solution of the differential equation (4). By the method of 
the preceding section, the characteristic polynomial z* + 100 has roots 
+i10, and the corresponding homogeneous differential equation x” + 100x 
= 0 has the general solution 


A cos (10t) + Bsin (102). 


According to Theorem 1, the general solution of the original differential 
equation (4) is the sum of the particular solution of the original equation 
and the general solution of the homogeneous equation. Thus the general 
solution of equation (4) is 


y(t) = Acos (10t) + Bsin (10t) — 9.8. 
Use the initial value y(0) = 2 to find A. 
2 = Acos(0) + Bsin(0) —-98 = A — 9.8, A= 118. 


Now compute y’(t), and substitute the given initial value y’(0) = 50 to 
find B. 


y@ = —10A sin (10t) + 10B cos (100). 
50 = —10A sin (0) + 10B cos (0) = 10B, B= 5.0, 
The required particular solution is thus 
y@) = 11.8 cos (10t) + 5.0 sin (10t) — 9.8, 
shown in Figure 14.7.2. 


In the remaining examples we shall concentrate on the first part of the 


problem, finding some particular solution of the given differential equation. In 
each case we could then find the general solution by solving the corresponding 
homogeneous equation and applying Theorem 1 as we did in Example 1. 


EXAMPLE 2 Find a particular solution of the differential equation 


y" — y' — 6y = 5 + 1827. 
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(0,2) 


Figure 14.7.2 Example 1 


Since f(t) is a polynomial of degree two, we guess that some particular 
solution y(t) is a polynomial of degree two, 


y(t) = K+ Lt + Me’. 
First we find the first and second derivatives of y(z). 
y =L+2Mt, y” = 2M. 
Next we substitute these derivatives into the given differential equation. 
y” — y' — 6y = 2M — (L + 2Mt) — 6(K + Lt + Mt’) 
= (2M — L — 6K) + (—2M — 6L)t — 6M?’ 
= 5+ 1827. 


In the last equation the coefficients for each power of t must be equal. There 
are three equations, for units, t, and ¢?. 


units: 2M —L—6K=5 
t: —2M —-6L=0 
t?: —6M = 18. 
We can now solve the three equations for the three unknowns K, L, and M. 
K = -2, E=1, M = -3. 
The required particular solution is then 


y(t) = —2 +t — 32’. 


It can be shown that whenever the forcing term f(t) is a polynomial of 
degree n, the differential equation (1) will have a particular solution that is a poly- 
nomial of degree n. When f(t) is a polynomial of degree n, the guess y(t) should be a 
polynomial of degree n with unknown coefficients, 


y(t) = Ag + Ait +--+ +A," 
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EXAMPLE 3 Find a particular solution of the differential equation 
y+ Ty + 10y = e%4, 
We guess that there is a particular solution that is a constant times e7 
y(t) = Me. 
The first two derivatives of y(t) are 
y(t) = 3Me*, y"t = 9Me*!, 

Substitute these into the original differential equation. 

9Me™ + 21Me* + 10Me* = e™. 


Cancel the e*', and solve for the unknown constant M. 
1 
9M + 21M + 10M = 1, M= io 0.025. 


The required particular solution is 


y(t) = 0.025e*, 


Here is the rule for guessing a particular solution of the differential equation 
of the form (1) when the forcing term f(t) is an exponential function f(t) = e“'. We 
first should find the roots of the characteristic polynomial az? + bz + c. If k is not 
a root of the characteristic polynomial, there is a particular solution of the form 
y(t) = Me™ (as in Example 3 above). If k is a single root of the characteristic poly- 
nomial, there is a particular solution of the form y(t) = Mte™. If k is a double root of 
the characteristic polynomial, there is a particular solution of the form y(t) = Mt*e™. 


EXAMPLE 4 Find a particular solution of the differential equation 
y" ae Ty" 4 10y = e7 2! 


The characteristic polynomial] z? + 7z + 10 has roots —2 and —5. Since 
—2 isa single root of the characteristic polynomial, our guess at a particular 
solution should be 


y(t) = Mte~*", 
The first two derivatives of y(t) are 
y(t) = Me~?* — 2Mte~*", 
y(t) = —4Me™*' + 4Mte~*#. 
Now substitute into the original differential equation. 
e7t = yp” + Ty’ + 10y 
= —4Me~*' + 4Mte~*" + 7Me~*' — 14Mte~*' + 10Mte~*! 
= 3Me*. 


Then M = 1/3, and the required particular solution is 


1 
y(t) = (F)ren*s 
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In this example the simpler guess Le~ ** would not have worked. The trouble 
is that Le” *" is a solution of the corresponding homogeneous equation, so it cannot 
also be a solution of the original differential equation. To see what happens, let us 
try to use the method with the guess u(t) = Le~ 4. Computing the first two derivatives 
and substituting, we get 


u(t) = —2Le-, u(t) = 4Le™, 
eo =u" + Ju’ + 10u = 4Le~?* — 14Le7 2" + 10Le7 2", 


The right side of the above equation adds up to zero, so we cannot solve for the 
unknown constant L. 

In a physical system, the forcing term is often a simple oscillation that can 
be represented by a function of the form f(t) = Gcos (wt) + H sin (wt). Here is the 
rule for guessing a particular solution of the differential equation (1) when the 
forcing term is f(t) = Gcos (wt) + H sin (wt). If z = iw is not a root of the charac- 
teristic polynomial, then the differential equation will have a particular solution of the 
form y(t) = K cos (wt) + Lsin(wt). On the other hand, if z = iw is a root of the 
characteristic polynomial, then there is a particular solution of the form y(t) = Kt cos(wt) 
+ Lt sin (wt). 


EXAMPLE 5 Find a particular solution of the differential equation 
yp’ + l6y = —sin (42). 


The characteristic polynomial has roots +i4. Then cos (4t) and sin (4t) are 
already solutions of the homogeneous equation, so our guess must have an 
extra factor of f. The guess for a solution is then 


y(t) = Kt cos (4t) + Lt sin (41). 
Compute the first two derivatives of y(t). 
y(t) = K[—4t sin (4t) + cos (42)] + L[4t cos (4t) + sin (42)], 
y"(t) = K[—16t cos (4t) —.8 sin (4t)] + L[—16t sin (41) + 8 cos (4t)]. 
Now substitute into the original differential equation. 
K[(—16t + 161) cos (4t) — 8 sin (4t)] 
+ L{(—16t + 162) sin (4t) + 8 cos (4t)] = —sin (42). 
—8K sin (4t) + 8L cos (4t) = ~Sin (42). 


From the sine terms we get —-8K = —1,so K = ¢. From the cosine terms we 
get 8L = 0, so L = 0. The particular solution is therefore 


y(t) = 0.125t cos (41). 


In Example 5, the particular solution oscillates more and more wildly as 
t approaches infinity, as shown in Figure 14.7.3. This happens because the forcing 
term cos (4t) has the same frequency as the solutions of the homogencous equation, 
Acos(4t) + B sin (4t). In this case the forcing term causes the oscillation to build up. 
This phenomenon is called resonance. 

If, instead, the forcing term in Example 5 had a different frequency, — sin(q?) 
where @ is not equal to 4, then the particular solution of the differential equation 
would be a simple oscillation of the form K cos (mt) + L sin (wt), whose amplitude 
does not change with time. 
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V(t) = 0.125 ¢ cos(47) 


Figure 14.7.3 Exercise 5 


EXAMPLE 6 Find a particular solution of the differential equation 
y" + 6y’ + 25y = cos (46). 
The characteristic polynomial z? + 6z + 25 has roots —3 + i4. Since i4 
itself is not a root, the proper guess is 
y(t) = K cos (4r) + L sin (42). 
Both a cosine term and a sine term are required, because the derivative of 
a'sine is a cosine. Compute the first two derivatives of y(t). 
y'@ = —4K sin (4t) + 4L cos (41), 
y(t) = —16K cos (4t) — 16L sin (42). 
Now substitute into the original differential equation. 
cos (4t) = [—16K cos (4t) — 16L sin (4t)] + 6[—4K sin (41) + 4L cos (41)] 
+ 25(K cos (4t) + Lsin (41)] 
= (—16K + 24L + 25K) cos (4t) + (—16L — 24K + 25L)sin (40). 
Both the sin (4t) coefficients and the cos (4t) coefficients must be the same 


on both sides of the equation. Thus we have two equations in the two un- 
knowns K and L. 


sin (4f) terms: OQ = —24K +4 9L. 
cos (4t) terms: 1=9K + 24L. 


Solve for K and L. 


I 8 
K=7 = Fig" 


The required particular solution is 


y= (3) cos (4t) + Gs) sin (41). 
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In Example 6, the particular solution is a simple oscillation, while the 
general solution of the corresponding homogeneous equation is a damped oscillation. 
The general solution is their sum: 


219 


The first term, e~ **[A cos (41) + B sin (4t)], approaches zero as t > 00 and is called 
the transient part of the solution. The second term, [(1/73) cos (4t) + (8/219) sin (41)], 
is called the steady part of the solution. The constants A and B, which depend on 
the initial conditions, appear only in the transient part and not in the steady state 
part of the solution. No matter what the initial conditions are, every particular 
solution will approach the steady state part of the solution as t > oo. The effect of 
the initial conditions dies out as f > ©. 

The same thing happens in any mass-spring system with damping where 
the forcing term is a simple oscillation. Let us consider a mass-spring system 


y(t) = e *[A cos (41) + Bsin (40)] + l(a) cos (4t) + Gs) sin «| : 


(5) my” + by’ + ky = cos (at), 


where m, b, and k are positive and the forcing term cos (wt) has frequency w. A 
particular solution can be found of the form 


K cos (wt) + Lsin (ot). 


The constants K and L can be computed as in Example 6. The general solution of 
equation (5) is 


y(t) = e “[A cos (Bt) + Bsin (Bt)] + LK cos (wt) + Lsin (at). 


As in Example 6, the first term approaches 0 as t > 00 and is called the transient part 
of the solution, and the second term is called the steady state part of the solution. 
Again, every particular solution of the mass-spring system will approach the steady 
state part of the solution as t > oo. 


PROBLEMS FOR SECTION 14.7 


In Problems 1-12, find a particular solution of the given differential equation. 


1 y” — 10y’ + 25y = cost 2 y" + 2./2y' + 2y = 10 
3 y” + 16y = 82? + 3t —4 4 y” + 2y = cos (5t) + sin (St) 
5 y” — 2y' + 2p =e” 6 y’ —6y+13=142r+e! 
7 yh ty —b6y=e™ 8 oy’ ay + By =e™ 
9 y” + 16y = cos (4t) 10 y” + 9y = 3 sin (30) 

11 y” + 12y’ + 36y = 6e°* 12 y” — 8y' + 16y = —2e* 

In Problems 13-16, find the general solution of the given differential equation. 

13 y" + 6y' + Sy = 4 14 y’—y—12p=t 

15 y’ + 5y = 8sin (22) 16 y” — 4y = 4e”" 


In Problems 17-20, find the particular solution of the initial value problem. 
17 y-y=3tts5, y(O) = 0, y'(0) =0 
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y+9y=4,  yOQ=0, yO) =0 

y" + 12y' + 37y = 10e7*, = KO) = 4, "(0) = 0 

y" + 6)" + 18) = cost — sine, y(0) = 0, y(0) = 2 

A mass-spring system mx” + bx’ + kx = 689 cos (2f) has an external force of 689 cos (25) 


dynes, spring constant k = 29, damping constant b = 4, and mass m = I gm. Find the 
general solution for the motion of the spring and the steady state part of the solution. 


A mass-spring system mx” + bx’ + kx = 2sint has an external force of 2 sint dynes, 
spring constant k = 24, damping constant b = 12, and mass m =3gm. Find the 
general solution for the motion of the spring and the steady state part of the solution. 


In the mass-spring system 
(5) my” + by’ + ky = cos (at), 


where m, b, and k are positive, show that the steady state part of the solution has 
amplitude 
1 


/(k — ma?)? + b2 
In Problem 23, show that the frequency w in the forcing term for which the steady 
state has the largest amplitude is 
jk 
arc, d m 


and the largest amplitude is 1/b. This frequency w is called the resonant frequency. 


JK? +1? = 


Using Problem 24, find the resonant frequency for the mass-spring system 
y" + 6y’ + 25p = cos (tr). 


PROBLEMS FOR CHAPTER 14 


Find the general solution of y’ = 3? cos t. 
Find the general solution of y" = ,/y. 
Solve y" = e*/t, x1) = 2. 

Solve y’ = P/Q" + 1), (0) = 1. 

Find the general solution of )” + 10ry = 0. 
Find the general solution of y’ + e“'y = 0. 
Solve y’ + 6y = 0, (0) = 1. 

Solve y’ — y\/t = 0, (0) = 2. 

Find the general solution of )’ + 3y = 2r. 
Find the general solution of )’ — ty = r?. 


A population has a net birthrate of 2% per year and a constant net immigration rate 
of 50,000 per year. At time t = 0, the population is one million. Find the population y 
as a function of t. 

Repeat Problem 11 for a net immigration rate of — 50,000 per year (that is, emigration 
exceeds immigration by 50,000 per year). 

Show that the initial value problem y’ = cos (y” + 8), y(0) = 1, has a unique solution 
forO<t< <x. 

Show that the initia] value problem y’ = 1/(2 + sin ), (0) = I, has a unique solution 
forO<t< <x. 


Find the general solution of y” — Sy" + 4y = 0. 


EXTRA PROBLEMS FOR CHAPTER 14 


Find the general solution of y” + 400) = 0. 

Find the general solution of y” — 4y’ + 8y = 0. 

Find the general solution of y” — 14y’ + 49y = 0. 

Solve y” + 4y’ — Sy = 0, (0) = 0, y'(O) = 1. 

Solve y” — 20y’ + 100y = 0, y(0) = 1, y’'(0) = 0. 

A mass-spring system mx” + bx’ + kx = 0 has mass m = 2 gm and constants b=6 
and k = 5. At time t = 0, its position is x(0) = 10 and its velocity is x((0) = 0. Find 
its position x as a function of t. 

Work Problem 21 if the system is subjected to constant external force of 3 dynes. 

Find the general solution of y’ — Sy’ + 4y =2 +1. 

Find the general solution of y“ + 400y = e’. 

Find the general solution of y” — 4y’ + 8y = cost. 

Find the general solution of y” — 14y’ + 49y = i: 

Solve y” + 4y’ — Sy = 26sin t, y(0) = 0, y'(0) = 0. 

Solve y” — 20y’ + 100y = e'™, y(0) = 0, y(0) = 0. 


901 


